Mapping trapped atomic gas with spin-orbit coupling to quantum Rabi-like model 
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We construct a connection of the ultracold atomic system in a harmonic trap with Raman-induced 
spin-orbit coupling to the quantum Rabi-like model. By mapping the trapped atomic system to a 
Rabi-like model, we can get the exact solution of the Rabi-like model following the methods to solve 
the quantum Rabi model. The existence of such a mapping implies that we can study the basic 
model in quantum optics by using trapped atomic gases with spin-orbit coupling. 
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I. INTRODUCTION 

Recently, many attentions have been paid to the cold 
atomic system with synthetic gauge field and spin-orbit 
(SO) coupling, which have been successfully realized in 
ultracold Bose gases [H-Q and Fermi gases @, Q. The 
experimental progress stimulated the intensive theoret- 
ical studies, including exploring schemes to create gen- 
eral gauge fields [l(| El and studying various interesting 
phases in these novel atomic systems with SO coupling 
[12| - tla |. While many theoretical works focused on the 
uniform systems with isotropic Rashba-type SO coupling, 
trapped systems with an external harmonic trap have 
also been addressed [l6l - [20| . It was shown that the pres- 
ence of a confining potential may qualitatively change the 
physical properties of atomic gases with SO coupling. 

In this work, we shall consider the anisotropic one- 
dimensional (ID) SO coupling realized in current exper- 
iments and scrutinize the problem of cold atomic system 
with SO coupling in a harmonic trap. We find that even 
the single particle problem of the trapped atomic sys- 
tem is highly nontrivial and not easy to be solved ana- 
lytically when both the SO coupling term and Raman- 
coupling term exist. As we shall show in the context, the 
trapped atomic system with Raman-induced spin-orbit 
coupling is equivalent to a quantum Rabi-like model. As 
a paradigm for modeling the simplest light-atom interact- 
ing quantum system [2lL |22| , the quantum Rabi model 
2l| was surprisingly not able to be analytically solved 
and only very recently has an analytical solution been 
found by Braak [23[. The theoretical progress has re- 
newed the interest in the study of the quantum Rabi 
model [2414281 . Particularly, Braak's solution can be also 
re-derived [24 1 within the extended coherent stats p9ll30j 
in a more straightforward way. 

Within the same scheme of solving the quantum Rabi 
model [HI, Hil , we can solve the Rabi-like model and get 
an analytical solution in terms of Braak's transcendental 
functions. By using the analytical solution, we can cal- 
culate the energy spectrum and dynamics of the Rabi- 



like model exactly. As the Rabi model has been widely 
applied to different fields of physics, including quantum 
optics, the cavity and circuit quantum electrodynamics 
[3TI HH and semiconductor systems [H], et al., we ex- 
pect that the mapping of the trapped atomic system with 
spin-orbit coupling to quantum Rabi-like model paves the 
way to study physical phenomena related to quantum 
Rabi model by using cold atomic systems. 



II. THE TRAPPED ATOMIC GAS WITH SO 
COUPLING 

The Hamiltonian of a single atom with SO coupling 
takes form 



H 3d = p 2 /(2m) + V so + V(f ), 



(1) 



with V(r) = ^m(uj 2 x 2 + w 2 y 2 + uj 2 z z 2 ) the trap potential, 
m the atomic mass, and p the atomic momentum oper- 
ator. The term of Vso describes the SO coupling. For 
effective spin- 1/2 systems considered in Refs. 0, [j], one 
has 



Vso = 2k r p x a z 



(2) 



where k r is the recoil momentum, f2 is the Raman- 
coupling strength, 6 is the two-photon detuning and a 
represents the Pauli operators. For the SO coupling de- 
scribed by Eq. ([2]) , the Hamiltonian can be separated into 
Hzd = H(x) + H2d(y,z) with the effective ID Hamilto- 
nian given by 



H = 



1 22 
H — muj x 



Vso, 



(3) 



2m 

where we have used uj — uj x for brevity. The model of 
([3]) has a deceptively simple form. However, as we shall 
show in the next calculation, the model is not easy to 
be solved except the special cases with either to = or 
ft = 0. 



III. MAPPING TO THE RABI-LIKE MODEL 
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To make a connection of the SOC model in the har- 
monic trap to the well known quantum Rabi model, we 
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use the representation of ladder operators, i.e., a and its 
adjoint a', to rewrite the Eq.Q as 

H = fojj(a)a + -) + 5a z /2 + VLa x /2 + ig(a) - a)a z , (4) 



where g = k r \/2mhu} and p x = iy / mhuj/2(a^ — a) is used. 
To simplify the above equation and compare with Rabi 
model, we make a shift H —> H — hu/2 and take h = 1. 
Now the Hamiltonian can be rewritten as 

H = uxJa + -a z + — a x + ig(a) - a)a z . (5) 

We note that a generalized quantum Rabi model can be 
represented as 

Hr = LocJa + ^a z + ^a x + g(a) + a)a z , (6) 

with 6 = corresponding to the quantum Rabi model. 
Comparing Eq.([5|) with the generalized Rabi model, we 
find that the only difference is the spin-boson coupling 
term. While the coupling term in the generalized Rabi 
model describes a spin-space coupling, the coupling term 
in the Rabi-like model of the trapped atomic gas de- 
scribes a spin-momentum coupling. 

The resemblance of these two models suggests us that 
the Rabi-like model ([5]) can be solved by applying sim- 
ilar methods to solve the quantum Rabi mod el p3l . |24| . 
Here we shall use the method of Chen et. al |2J| to de- 
rive our solution of Rabi-like model. For brevity, we give 
results and key steps of derivation in the main text, but 
leave some details of derivation in the appendix. To di- 
agonalize the Rabi-like model, it is instructive to span 
the Hamiltonian ([5]) in the spin space with the following 
matrix form 



H = 



a^a + ig (a 

Q 
2 



t _ 



n 

2 

a^a — ig (a 



t _ n\ - £ 



n (7) 

Here, in order to keep consistent with Ref. [24]], we take 
the matrix form in units of u> = 1. By introducing the 
Bogoliubov transformations A = a + ig and B = a — ig, 
the diagonal matrix element Hu and H22 can be diago- 
nalized as A' A and B^B, respectively. 

In terms of operator A, we can rewrite the Hamiltonian 

as 



H 



where 



A^A- 
2 



n 
2 



A*A - 2ig (At -A) +/3 



(8) 



2 S 9 5 

u = g - -, p = ig" - -. 

To diagonalize the Hamiltonian, we take the wavefunc- 
tion as 



|V.) = J2 > E ^ J< n \n) A , (9) 



where \n) A = \0)a ^ s the extended coherent state 

with \0) A = e-2 9-- l ° a |0) o . The expansion coefficients 
J~ and K~ can be determined self-consistently by solv- 
ing the eigenvalue equation H \ip) = E\ip), which gives 
some restricted relations to coefficients J~ and K~ [24| 
(see appendix). 

Similarly, in terms of operator £?, the Hamiltonian is 
written as 



H = 
where 



= ( B^B + 2ig (fit _ B) + p § 

B^B - a' 



2 



(10) 



The wavefunction can be also represented in terms of B 
as 



,£(-l) n V^J+|n) B 

n=0 / 
(11) 



(fit)" 

where \n) B = ^-j |0) B is the extended coherent state 

with |0) B = e-y 2 ^ |0) a , and K+ and J+ are the 
coefficients to be determined by solving the eigenvalue 
equation H \t/j') = E \ip'). 

Since both wavefunctions \ip) and \ip') are eigenfunc- 
tions for the same eigenvalue E, they should be different 
by a complex constant if they are not degenerate. This 
requirement leads to a self-consistent condition between 
the two sets of coefficients, i.e., 

00 00 00 00 

E J n(-wT E J n Hs) n = E K ni-ig) n E *CHs)", 



n=0 



n=0 



71=0 



n=0 



The above equation can be rewritten as a compact form 
Gs(x) = with the transcendental function defined as 

Gs(x) = (^) 2 R + (x)R~(x) - R+(x)R-(x), (12) 
where x = E + g 2 

00 

r ± (x) = J2 K tW HaT, 

71 = 



R^x) = E 



77 = X - 71 ± 2 



(-ig) n - 



The K^{x) are defined recursively, 



nKt = fti(x)Kti+Kt2 



(13) 



with the initial condition K, 



- 1, Kf = f±(x), and 
(fi/2) 2 



n-x±6/2 



4(i-fiT V 2 ) 



(14) 



3 



,x 1(T 



u 



-1 01234567 




FIG. 1: (Color online) The energy of the Rabi-like model can 
be determined from the points of intersection of G-function. 
Here we set Q = 1.4, g = 0.1, 5 = and made a truncation of 
first 199 coherent state orbits. 

The eigcncncrgies can be determined by the zeros of the 
transcendental function Gs(x) [23|. For S — 0, Gg(x) = 
can be reduced to G (x) = with 

Gt(x) = J2 K n (x) ( 1 ± ) (-ig) n . (15) 



IV. RESULTS AND DISCUSSIONS 

The spectrum of the Rabi-like model can be deter- 
mined from the zeros of transcendental function (fT2|) or 
(fl"5j) . As illustrated in Fig.l, we can read the eigenenergy 
from the intersection points of G-function with horizontal 
axis (xi n ) via E — Xi n — g 2 . Using the above G-function, 
it is quite easy to get the relation of energy spectrum 
as function of coupling constant. To give a concrete ex- 
ample, we plot the energy spectrum for the system with 
SI = 1.4 in Fig. 2. As shown in the figure, the ground 
state energy decreases with increasing coupling constant 
g. When 6/2 is not a multiple of ui/2, no level crossings 
appear. This is consistent with Braak's results (23j . 

We note that the G function contains all the informa- 
tion of our system, including both eigenvalues and eigen- 
functions. As the eigenvalues are determined, coefficients 
of the corresponding eigenfunctions can be obtained by 
recursion relations (|13|) . Although the G-function is rep- 
resented as a summation of infinite power series, un- 
doubtedly practical calculations of the zeros of G function 
require truncations of coherent state orbits. Neverthe- 
less, the compact form of our solution implies that we can 
get desired accurate results in our scheme. A connection 
with the quantum Rabi model inspires us to investigate 
the dynamical behavior of the Rabi-like model, which 
may stimulate the study of quantum dynamics related to 
quantum Rabi model by using cold atomic systems. 

In order to study the dynamics of Rabi-like model, we 
first prepare the bosonic atoms in a harmonic trap and 



FIG. 2: (Color online) The energy of the Rabi-like model 
versus coupling strength g. Here we set Q = 1.4, 5 = 1 for (a) 
and 6 = 1.4 for (b). In the calculation, ui = 1 has been taken 
as units. 

suppose that the interaction between atoms is very weak 
and the system can be treated as noninteracting bosons. 
In principle, one can tune the interaction between atoms 
to be zero by the Feshbach resonance. We take the initial 
state as |0) o <£> (1,0) T , which is the ground state of the 
system in a harmonic trap. At time T = 0, we turn 
on the Raman-induced SO coupling. The system will 
evolve under the Hamiltonian ([5]) and the time-dependent 
wavefunction reads 

n \Wn \iPn) 

where \ip n ) is the eigenstate of the Rabi-like model cor- 
responding to the eigenvalue E n . The analytical solution 
of the Rabi-like model enables us to calculate the time 
evolution of the expectation value of <j z in a numerically 
exact way [HI, HH . 

The expectation value of a z is defined as P z {t) = 
(ip(t)\a z \ip(t)}, which describes the population differ- 
ence of two component atoms. The left panels of Fig. 3 
show the time-dependent population difference P z (t) for 
SI = 0.5, 5 = under various coupling strengths. For 
the weak coupling strength with g = 0.1, P z (t) oscillates 
between 1 and — 1, which indicates the periodic popula- 
tion swapping between different component atoms. With 
increasing the coupling strength, the typical Rabi oscil- 
lation in the weak coupling regime breaks down. For 
g = 1.5, P z {t) displays a sharp steplike decay with sim- 
ilar dynamic behavior as in the quantum Rabi model 
[28j . In the right panels of Fig. 3, we show the time- 
dependent population difference for system with SI = 0.5 
and 5 = 0.4. It is obvious that the oscillation amplitude 
of P z (t) is suppressed by the nonzero term of 5. From 
the Hamiltonian ([S]), given the initial state as (1,0) T in 
a z representation, one may understand the competition 
between the SI term and 5 term. While the nondiagonal 
term of SI flips the spin, the diagonal term of 5 tends 
to stabilize the initial state. The effect of coupling g is 
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somewhat like the diagonal term while in a more subtle 
way. 

In order to compare with the well-known quantum 
Rabi model, we take f2 with the same energy order as 
w, which is generally fulfilled for cases of quantum opti- 
cal systems. However, the energy scale of cold atomic 
systems is quite different from quantum optical sys- 
tems. In quantum optical systems, !1 « w and the 
coupling g is the smallest quantity, while in cold atom 
systems, the harmonic trap is usually not deep enough 
and generally Q 3> uj. For example, in the recent ex- 
periment of USTC group @, one has U) = (2ir)50Hz, 
Erecoii = (2ir)2.21kHz, S = AE recoi i, Q sa E reco u, from 
which we can get f2 sa 40 and 8 s=s 160 in units of u>. To 
see the effect of a large il, in Fig. 4, we give an example 
for the case with il = 10. Comparing with Fig. 3, we 
can clearly see the oscillation frequency increases with 
increasing the Raman coupling strength. The oscillation 
amplitude is strongly suppressed by the d term. 

The mapping of the trapped atomic gas with spin-orbit 
coupling to quantum Rabi-likc model may stimulate the 
study of interesting phenomena related to the basic quan- 
tum optical model by using the cold atomic system. Al- 
though the parameters of trapped cold atomic systems 
in current experiment are much different from the gen- 
eral Rabi model, the good tunability of cold atomic sys- 
tems may make it be possible to access the parameter 
regime of traditional quantum Rabi model by increas- 
ing the frequency of trap potential. One possible way is 
using the ID deep optical lattice to produce a series of 
two-dimensional (2D) Bose gases. In the direction of ID 
optical lattice, each 2D gas can be viewed to be tightly 
trapped by an effective harmonic trap. By this way, the 
frequency of the effective harmonic trap may be tuned 
to the same order of Raman-coupling strength [3|| . As 
quantum Rabi model in the quantum optical system and 
SO coupled atomic gases can be viewed as two limit cases 
in parameter spaces, one may also explore some novel 
properties of the Rabi-like model by using cold atomic 
systems in the parameter regime which is not accessible 
in traditional optical systems. 



SUMMARY 




FIG. 3: (Color online) Expectation value of a z for 5 = (red 
line) and S — 0.4 (blue line). From (a) to (c) or (d) to (f), 
the coupling strength is taken as 0.1, 0.7 and 1.5, respectively. 
Here $1 = 0.5, time t is in units of 2n/u) 
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FIG. 4: (Color online) Expectation value of a z for 5 = (red 
line) and 5 = 40 (blue line). From (a) to (c) or (d) to (f), the 
coupling strength is taken as 0.1, 0.7 and 1.5, respectively. 
Here J! = 10, time t is in units of 2n/u 



In conclusion, we show that the ultracold atomic sys- 
tem in a harmonic trap with Raman-induced SO coupling 
is equivalent to a Rabi-like model and the closed form of 
the exact solution is obtained. The connection of the 
trapped cold atoms with SO coupling to the Rabi-likc 
model enables us to explore a variety of properties of 
the cold atomic system by using traditional methods in 
quantum optics. We want to point out that this connec- 
tion is rather implicit as the harmonic trap here provide 
a source of bosonic field, which was only considered as 
a confinement in previous studies. Another reason of 
this implicity is the different energy scales in quantum 
optics and cold atoms. Rapid progress in cold atomic ex- 



periments gives the possibility to prepare our system in 
quantum optics limit. Our study also provides a new vi- 
sion for understanding the trapped cold atomic systems 
with SO coupling. 
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Appendix A: Derivation of solution. 

For simplicity, we shall use the method of extended 
coherent states developed by Chen et. al. (24[, which 
was used to re-derive the Braak's solution to the quantum 
Rabi model. Following the procedure for the Rabi model 
[IH , we can get the solution of our Rabi- like model. 

Left-multiplying a (n\ to the eigenvalue equation 
H\ip) = E\tf>), we can get the following restricted re- 
lations to the coefficients J~ and K~: 



(n-a-E)J n 



(n + [3~E) K- + 2ig (n + 1) R- +1 - 2igK~_ 1 
- 2 J n ■ 

From the above equations, we can determine the coeffi- 
cient f n recursively, 



with K+ = 1 and K+ = /+. 

As alternative representations of a non-degenerate 
state with the eigenvalue E, the wavefunction ([9]) and 
(JTTJ) should be equivalent, i.e., \ip) = r \ip'), where r is a 
constant. The above equivalent requirement leads to the 
following relations: 



oo 

E 

n=0 



Vn\n) A = r^(-l) n VriK+\n) Bi (A7) 



\K~ \n) A = r ^(-l)VnL# \n) B . (A8) 

n=0 n=0 

Left multiplying the vacuum state a (0| to both sides of 
the above equations and making use of the relation of 

Vri. a (0\n) A = (-l)V^! o (0|n) B = e-«*/ 2 (-ig) n , we can 
eliminate the ratio constant r and get the following rela- 
tion 



Jn 
nK~ 



-K~ 



2{n-a-EY 



Jn 2ig 



(n + p-E)- 



2 



A(n-a-E) 



(Af) 
(A2) 
(A3) 



with Kq = 1 and = 

Similarly, proceeding as before, we get the relations for 
two coefficients A'+ and J+: 



2 " 



(n - a' - E) J+ 
(n + p'- E) K+ + 2ig (n + 1) K+ +1 - 2tgK+_ 1 



' J n 



then we have 
/+ - 

nK+ = 
r+ _ 



n 



r = — 

Jn 2ig 



2(n - a' -E) 

•4 

n-2> 

(7i + p-E)- 



f + K + 4- K + 
-1 



n 2 



4 (n - a' -E) 



(A4) 
(A5) 
(A6) 



71=0 n=0 
oo oo 

^^(-iffr^ir+f-ifl)-. (A9) 



n=0 



n=0 



Making substitutions in terms of Eqs. (|Af I) and (|A4|) . we 
obtain 



E 



fi/2 



n/2 



n-a-E ny t-^n-a'-E 

n=0 n=0 



J2K-(-ig) n J2K+(-igy 

n=0 n=0 



(A10) 



In terms of the the transcendental function G$(x) de- 
fined by Eq. (|f 2[) in the main text, the above equation is 
nothing else but Gs(x) = with x = E + g 2 . 
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